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Abstract 

Superconductivity in small metallic specimens is studied with regard to the size 
dependence of the parity gap (Ap), a parameter distinguishing between the energy 
of even and odd number of electrons in the granule. Ap is shown to be an increasing 
function of level spacing 5. The energy gap of superconductor A, on the other hand, 
decreases with increasing S and vanishes at 5 = 8 C which is of the order of A. How- 
ever, nonzero value of Ap persists above 5 C in a gapless superconducting- insulating 
state. Level degeneracy in small specimens having perfect geometry changes the size 
dependence of the parity gap, the Josephson effect, and flux quantization. Parity gap 
is evaluated using an interpolation procedure between the continuum limit < A), 
the moderate mesoscopic regime (5 ~ A), and the nanoscopic scale (5 S> A), for 
which an exact solution to the pairing problem is provided with the numeric diag- 
onalization of system Hamiltonian in a small metallic cluster. 
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1 Introduction 



Superconductivity in metals is related to condensation of electrons into pairs 
which are extended objects with a size of order £ ~ Uvf/Aq. A fundamen- 
tal question [1] is whether pairing can survive in small specimens with size d 
less than £. In energy scale this corresponds to condition that 5i > A where 
8\ = hvp/d is typical scale of dimensional energy quantization. The energy 
levels can remain highly degenerate in a perfect sample (see Fig. 1) but in 
a disordered sample, level spacing <5 2 ~ Ep/N s , where N s is the total num- 
ber of sites (which is of the order of total number of electrons, N), is much 
smaller than 5i and can put a limit to superconducting pairing at 5 2 ~ A . 
In recent experiments[2,3], signature of superconducting pairing was traced in 
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Fig. 1. Degeneracy of energy levels for parabolic energy vs momentum dependence 
of total energy E n = Cn 2 = C(n\ + + n|). 

ultrasmall specimens down to sizes of 10 nm. The superconductivity reveals 
itself through the number-parity effect, the dependence of superconducting 
energy on whether the number of electrons in sample is even or odd. Matveev 
and Larkin [4] introduced a parity gap Ap = ^(E 2n +i + E 2n -i) — E 2n where 
En is superconducting energy vs number of electrons N, to characterize that 
the paired state is preferable over the state with odd number of electrons. In 
a number of papers, in particular in the paper of Braun and Von Delft [5], 
E N was calculated assuming that in an odd-number sample, electron state 
near the Fermi energy is blocked from pairing, similar to the way in which 
finite-number systems are treated in nuclear physics [6]. In present paper, we 
investigate the suppression of superconductivity in a finite size system by us- 
ing an approximation similar to that of [5], and also in the opposite limit, 
by starting from atomic-scale system (small atomic cluster) for which an ex- 
act solution of pairing problem is provided, and trying to match small- and 
large-size approximations through an interpolation procedure. 



2 Large-scale approximation 



In the Bardeen-Cooper-Schrieffer (BCS) approximation [7] a variational den- 
sity matrix p Q/ g = 4>0* exp(— (H — fiN) /kT)^^ is exploited to obtain the 
ground state energy with trial wave functions 

= ITK + ^k4 T al ki )|0> and = J] («k + v^al H )al F \0) (1) 



for even and odd number of electrons, respectively. In the BCS pairing scheme 
(case 1), assuming equally spaced levels with distance 8, we receive A 2m = 
^2m+i = A whereas in the blocked pairing state (case 2) A 2m = A and 
A 2m +i = A — 5/2. The energy of superconductor at T = is calculated as 
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Fig. 2. Josephson critical current vs energy level spacing 5 (solid line) and gap 
parameter A vs 5 (dashed line). The latter has been scaled by a factor of 1/2 for 
convenience. 

(H — fiN) = J2n{£,n — \J A 2 + £^) resulting in a parity gap value A^P = | (for 
5 <C A ) in the first case, and Ap" 1 ~ A In + |<5 in the second. Energy 
of the blocked state is larger than that of the BCS state and therefore there 
is no reason to believe that the former is a better approximation for pairing. 
There is no direct relation of the above calculation to the question that there 
is a large uncertainty in the number of particles, and assertion that number 
of particles in 0^ is odd is merely an illusion. 



3 Josephson effect in mesoscopic metal particles 

Consider two identical small superconducting grains interacting through the 
Hamiltonian H T = J2nmTnm a nb m + H.c. where at(a n ) creates (annihilates) 
an electron in the first grain in state n. Similarly and b n are operators for 
the second grain. By introducing the phases <f>i and 2 to the Bogoliubov u, v- 
coefficients in Eq.(l), we calculate coupling energy due to as a function of 
4>i and 02 hi a standard way [7] AE = AE' S — AE$ cos(0i — 2 ) where 

AE S = (\T\^ £ 1 ■ (2) 

\tn\,\Z m \«*D 1 tm > tntm 

The stationary Josephson current at zero temperature is J c sin(0x — <f> 2 ) where 
J c = 2eAE s /h. Figure 2 shows dependence of Josephson critical current and 
energy gap as a function of gap suppression parameter 5/A . The curves show 
that superconductivity is suppressed at level spacing 5 C = 3.56A . This how- 
ever is a conclusion within the mean field approximation. It does not exclude 
the possibility that pairing may survive for 5 > 5 C as a gapless superconduct- 
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Fig. 3. Parity gap vs level spacing for nondegenerate case (a) and for double de- 
generate levels (b). In (a), numbers 1, 2, 3, 4 label dependences received with a 
number of levels in the interaction shell np, =30, 60, 120, and 360. In (b), num- 
bers 1, 2, 3, 4 correspond to index i in various pairing gaps defined according to 

ing state or as a state of superconducting insulator with an effective insulator 
gap S larger than the superconducting pseudogap. 



4 Richardson- Sherman approximation 

A finite size solution to the pairing problem was suggested by Richardson and 
Sherman [8] (and recently revived by von Delft [9]) assuming that states with 
energy within the interaction shell — uj d < £ n < uj d can be approximated as 
hard-core bosons whereas states away from this shell are treated exactly as 
non-interacting particles. The interaction Hamiltonian for hard core bosons is 

V5 

H = 26j2 nC L C na- -T C L C m/3 (3) 

na \£ n \,\U\«OD,al3 

where boson operators are c na = a„^a na -|-, d is level degeneracy, and a labels 
states within the degenerate energy level £ n = nS. Hamiltonian 3 is solved [10] 
by iteration scheme proposed by Richardson and Sherman [8] . Dependence of 
Ap on S/Aq is shown in Fig. 3. In nondegenerate case (d = 1) Ap(S) shows 
a behavior similar to previous results with smaller number of levels [11,12]. In 
case of degenerate energy levels, pairing gap is multi-valued depending upon 
which three successive states are used in calculation of Ap. The true pairing 
gap is possibly the largest one in Fig. 3(b), and others represent a new kind 
of mesoscopic effect which may or may not be related to superconductivity. 
Small 5 of Ap(S) does not match with any of pairing gaps calculated in Section 
3, which means possibly that neither of the approximations (original BCS or 
BCS with one level blocked) are good for small systems. 
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Fig. 4. (a)Energy vs number of electrons in cubic cluster with attractive interaction 
at sites. (1) U/t = 2, (2) U/t = 3, (3) U/t = 4. (b)Parity gap vs level spacing in 
cubic cluster. 

5 Parity effect in atomic clusters 

We calculate the exact eigenstates of the negative-?/ Hubbard Hamiltonian 

N 

H = -t ]T a\ a a ja ~UJ2 alAi a U a ^ + J2 v i a la a ia (4) 

<i,j>,<7 i=l icr 

(with U > 0) in small system of iV sites. Vi is random potential at site i. The 
fact that Hamiltonian 4 commutes with the total number of electrons and 
also with the total spin leads to lowering of dimensionality 2 2N x 2 2N of the 
matrix. Interaction is non-retarded (unlike BCS Hamiltonian) and momentum 
dependence of energy is nonparabolic. We do not believe however that these 
are crucial differences, and indeed the small-f/ limit mean-field solution [13] 
corresponds to an energy gap Ao — 12t exp(— Vlt/U). Figure 4, a shows mini- 
mal energy of cubic cluster as a function of number of electrons N. Parity 
effect clearly manifests itself in lower values for even iV and, when plotting 
A P against 5/A (Fig.4,b), displays a behavior very similar to one received in 
the Richardson- Sherman model. 
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